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ABSTRACT 

A multiple magnet bearing can be developed from N individual electromagnets. This paper presents the 
constitutive relationships for a single magnet in such a bearing. Analytical expressions are developed for a magnet 
with poles arranged circumferencially. Maxwell’s field equations are used so the model easily includes the effects of 
induced eddy currents due to the rotation of the journal, Eddy currents must be included in any dynamic model 
because they are the only speed dependent parameter and may lead to a critical speed for the bearing. The model 
is applicable to bearings using attraction or repulsion. 

INTRODUCTION 

Research and development activity on passive, active, and combination magnetic bearing systems spans over 150 
years beginning with Earnshaw in 1839. The authors will not attempt to summarize the history for this workshop 
but will simply note research in magnetic bearings, magnetic dampers, and magnetic levitation for vehicles. High 
speed magnetic bearings are a commercial reality, being used in grinding and polishing machinery, vacuum pumps, 
compressors, turbines, generators, and centrifuges. 

Passive systems using permanent magnets in repulsion are usually incapable of generating sufficient load 
carrying capacities. Two electromagnetic levitation methods have met with success: direct, position feedback 
control techniques; and ac modulated or indirect feedback methods. The latter suffers from high eddy current losses 
and a small range of stable air gaps. Early experiments with fully active systems (1957) were largely thwarted 
by the then high cost and large size of control system components. Since active magnetic bearings provide both 
damping and support, and with the reduction in cost and size of sophisticated electronics, the choice is clear. 
Subsequent efforts for load bearing situations have concentrated on active magnetic bearings. 

However, the dynamics of direct feedback magnetic bearings have not been as thoroughly analyzed as have 
journal bearings. Most of the available literature deals with empirical ideas and concentrates on reliability of 
the bearing, reducing the size, weight and complexity of the devices. The authors have completed an extensive 
investigation of the dynamics of an active magnetic bearing [1,2, 3, 4] in terms of equilibrium points, transient 
response, on-set of instability, limit cycle size, and forced response. The only speed dependent effect was the change 
in forces due to eddy currents. This change leads to a critical speed, a Hopf bifurcation with subcritical unstable 
limit cycles. The reader is referred to [3,4] for a description of the dynamic analysis and results. 

The purpose of this paper is to present the eddy current analysis. Although a good analytical model has not 
previously been available for eddy currents due to shaft rotation, a number of authors have calculated the eddy 
current effects in other geometrical configurations using a hypothetical simplified model and finite element methods 
[5,6]. Studies on linear induction motors are available [7,8] which can be extended to magnetic bearings by several 
assumptions and manipulations. In the above studies, the current density vector is cross multiplied to calculate 
the induced forces, which does not take care of the attractive force when the moving material is ferromagnetic. In 
this paper, the forces are calculated using the Maxwell’s stress tensor approach which in one calculation gives all 
the forces involved. Matsumara [9] has derived the fundamental equations for a horizontal shaft magnetic bearing 
taking into account the rolling, pitching, and yawing of the rotor. In deriving the equations of motion, he assumes 
that the journal consists of a laminated core and consequently no eddy currents are generated in the material. 

FOUR MAGNET BEARING 

Consider a bearing using four magnets as shown in Fig. 1. The shaft has high permeability and high conduc- 
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tivity and is not laminated. (Lamination would serve to inhibit eddy currents.) Purely for simplicity, all of the 
four electromagnets are identical, modelled as a coil of N turns on a laminated core of high permeability. Each 
magnet applies a flux density Bi and the dead weight is suspended by the difference in forces exerted by top and 
bottom magnets. All magnets generate a flux density which can be adjusted to vary the equilibrium gap lengths 
under the magnets. This fully energized configuration was chosen so that the magnets could generate the effect of 
repulsion by decreasing the attraction. The shaft displacement is measured by two coordinates (£,* 7 ) as shown in 
Fig. 1, measured from the center of the clearance space in the horizontal and vertical directions respectively when 
the journal is spinning under the magnets. 

If the shaft is not laminated, motion of the conducting shaft through the supporting magnetic fields will induce 
eddy currents. These eddy currents create two kinds of force on the rotor: drag forces which lead to additional 
power dissipation and coupling of motions of the journal in two perpendicular directions; and repulsive forces which 
tend to counter balance the attractive forces and shift the equilibrium point. 

The following geometrical assumptions are made. The journal always remains perfectly aligned within the 
bearing (no tilting). Under small displacements the surfaces of the journal and the magnet pole faces are assumed 
to remain parallel. Since the individual poles are located at angles ±0 relative to the control of the magnet, it is 
assumed that when the rotor moves vertically a distance 77 , the change in gap length for the vertical magnets is 
rjcosO. Similarly, when the rotor moves horizontally by a distance £, the change in gap length for the horizontal 
magnets is £co$0. Any other translational motion of the rotor can be written as a superposition of the motions 
in £ and 77 directions. The effect of unequal gap lengths under a vertical magnet caused by a rotor motion in the 
horizontal direction or vice versa is neglected because the total gap length under that magnet remains constant. 

EQUATIONS OF MOTION 

The differential equations of motion describing the response of the rotor system require the vectorially combined 
magnetic forces from all the magnets at a station along the rotor axis. 

= /(£,»?>£>»?. B l ,B 2 ,B 3 ,B 4 ) (1) 


Most rotordynamics codes linearize this to provide the eight bearing coefficients. 


(F(\ = \ C aC(r,]((\ , 
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( 2 ) 


It is immediately obvious that the dynamic characteristics of the bearing are determined by the control system 
which produces the changes in magnetic flux density. This system may be simply local, represented by 


Bi = f( Ai,A f ) (3) 

or coupled 

Bi = (4) 

The design of the control system will not be discussed here. Both pole placement methods [4,10] and linear 
quadratic regulator theory [ 1 ] can be used to design state variable controllers. 

Each magnet force will be expressed as resolutes normal and tangential to the magnet center line. These forces 
will include steady state attraction forces and eddy current forces (both repulsion and drag). Only spin velocity is 
assumed to generate eddy currents. Motion of the journal center (£,*?) doesn’t generate eddy currents. The force 
resolutes are only functions of flux density, gap, and journal speed. 


F ni = f( 

(5) 

F ti = /(A,,n,50 

( 6 ) 

Each magnet gap can be determined from the journal location 

A , = /(£,»/) 

(7) 
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It can then be shown that unless the control is based on velocities 


— C £ r ) ~ — C rjrj — 0 ( 8 ) 

and the eddy current drag forces produce bearing cross stiffness terms. 

Assuming the controller is determined, the general procedure is to: 

(a) Determine the equilibrium location for a given spin speed (This assumes no journal velocity, £ = 77 = 

°) ; 

(b) Linearize magnetic forces about equilibrium locations; 

(c) Form system dynamic equations, extract eigenvalues, determine stability; 

(d) Perform parameter study to locate system bifurcations; 

(e) Apply center manifold theorem to determine large amplitude behavior of any bifurcation points; 

The usual parameter study would be to follow the equilibrium location and stability as functions of spin speed. 
Note that due to nonlinearities, the bearing coefficients must be determined at each different parameter value. 


MAGNETIC FORCES INCLUDING EDDY CURRENT EFFECTS 

The complete eddy current analysis for 4 magnets as shown in Fig. 1 with unequal gap lengths is analytically 
intractable. (It is probably amenable to finite element techniques.) The following approximations are required. 

First, the problem is assumed two dimensional. 

Second, the problem will be unwrapped and considered as periodic on a half-space. However, before net forces 
are calculated, the surface tractions predicted by Maxwell’s stress tensor will be wrapped around a circular shaft. 
As the shaft starts spinning, the eddy currents tend to repel the applied magnetic field and the skin depth of 
penetration becomes very small. This motivates a semi infinite assumption in the radial direction. The velocity 
will be treated as constant with respect to y, although it actually decreases linearly as the center is approached. 
V = fl(y -f A - y). This assumption is supported by the very small skin depth. 

Third, each magnet will be considered separately, and the magnetic field for each magnet will be determined 
individually. The net force of each magnet is then determined, leading to 4 vector forces which are then summed 
vectorially. An alternative (more complex) solution is discussed later in the paper. Fig. 2 shows the approximate 
problem to be solved for each magnet. 

Fourth, the square wave applied flux density is expanded as a Fourier series 


/ x r 2b i 2cii , . 2 ix 

»(*) = Is ) ~ cos (-^-)] sir, (-j9-) 


( 8 ) 


or 

B m {x) = JZ Bisin(kix) (9) 

» 

The field density and current density distributions within the moving material solve a linear problem, and hence the 
principle of superposition can be invoked and each harmonic handled separately. Figure 3 shows a single harmonic 
applied to the semi-infinite medium. 

The equation that describes the distribution of the magnetic field in the conducting medium is derived from 
Maxwells field equations. 


1 9 - dB - 

D+ ar = r> ‘< v ' xB > < 10 > 

where B = B x e x + B y e y (e x and e y being unit vectors in x and y directions respectively). The y component B y 
is determined from the y component of (10) The remaining component B x can be determined from the relation 
V • B = 0. The magnetic field is driven by the applied magnetic field density, and so solutions with the same 
traveling wave dependence on ( x,t ) are assumed. That is, il is assumed that the flux density takes the form 

B = [B x (y)e x + B y ( !i)i y }c }f - tr ~ w B ( 11 ) 

The solution form in the y-direction is e ±qy where 


7/ 


i — y/T+js — /.*, « /l -h j 


jia \ ' 


A, 


(12) 
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Eqn (12) can be used to predict the skin depth. 

The solution domain is divided into two regions, denoted by subscripts 1 and 2 respectively: Region (1), the 
air gap where <r = 0(0 < y < A); and Region (2), the moving conducting medium (A < y < oo). The solutions 
within each region have two constants of integration, determined from the following 4 boundary conditions: 

(1) the applied flux density at y — 0; 

(2) the solution cannot grow as y — > oo (one constant is zero); 

(3) conservation of flux at the interface y - A (V • B = 0, using the divergence theorem, this leads to 

By i = Byi), 

(4) at the interface, V x H — 0 (using Stokes theorem, this leads to H x \ — II £ 2 )- 
Hence, the flux density distribution throughout both regions can be determined. 

In the air gap, the solution is 


j B u = (( D ie k ' ,J + D 2 e~ k ' y )e x + (£>ie l ' s - D 2 t- k ' v )e y )e^ k ' x ~ wt) 

(13) 

i 

j where 


! _ 

; 1 2[cosh(kiA) + ‘^■^■sinh(kiA)] 

(14) 

j (1 + 

j ° 2 ~ 2\cosh{kiA ) + ^f-sinh{ki A)] 

(15) 

In the journal, the solution is 


B-n = (^-e r + e v ) D 3 e _, ' y e j(t,x_ “' <) 

(16) 

where 


B i e q '^ 

1 ° 3 ~ [cosh(ki A) + ^sinh(kiA)] 

(17) 

The forces acting on the conducting medium are calculated by Maxwell’s stress tensor, 
with currents and no charges, the forces acting on a body are given by 

For magnetic problems 


F= f -{UB n- ]-B 2 n}dA (18) 

Jn t 1 2 

where E is any closed surface surrounding the body and not containing any other body and B is the value of the 
field on the closed surface. Choosing a closed surface E such that it extends from — oo to oo and includes only the 
conducting medium, the integration is carried out with n = -e y . 

The complete flux density distribution in the whole region of Fig. 2 due to all the applied sinusoidal waves is 
determined by superposing the individual fields. Each component of the field ( B x , B y ) is an infinite series in sine or 
cosine terms. The value of the B field at the interface, which is required for calculating the forces, is calculated by 
substituting y = A. It is interesting to note that the integral and the summation are interchangable in order. That 
is, the force for each component of the Fourier series can be determined and then summed or the fields summed 
and the force determined. This is perhaps surprising because the problem is nonlinear, but the infinite senes for 
B x and B y are made up of sine and cosine terms which are orthogonal to one another and all the cross terms drop 
out during integration over one period. 

The total flux can also be calculated and compared with that predicted by magnetic circuit theory (Appendix 
A). For the numerical example which follows, the more detailed solution is about 8% lower, showing the effects 
of magnetic circuit assumptions (uniform field density and no leakage in air gap). For comparison to classical 
problems, the forces for the problem in Fig. 2_will be determined, although the reader is cautioned to go on Eqn. 
(23) for a bearing problem. Substituting for B and evaluating tin integrals, the drag and lift forces per unit area 
acting on the material turn out to be 

p. = Liny , J Jl— * (i9) 

ras 4 /Iq t l I ( n ii + n 2 i ) Be\/l + 
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where 


4^0 V ("?< + "!<) v /* , v“ r ' <J 


and ni and n 2 are the real and imaginary parts of ji. 


ki 


K™) 


( 21 ) 


7i = [eosftffc.A) + — |i-sinA(jbiA)] (22) 

/i A: t 

As expected, when no currents are induced in the slab (V=0,s=0), there is no drag force and the lift is the 
magnetic attractive force. There is an optimum value of V at which the maximum force per unit area F x is produced. 
The lift force decreases as V is increased, at some value of V the lift force becomes zero, and at higher values the 
lift force acts in the opposite direction (repulsion). At very high values of V there is no drag because all the flux is 
excluded from the material and the repulsion force reaches an asymptotic value irrespective of the permeability of 
the material. 


For the magnetic bearing, the forces acting on a rotating shaft are calculated by wrapping one period of the B 
field distribution back onto the circular shaft. Choosing a closed surface £ on the circumference of the rotor and 
simplifying the integral in Eq. (18) give the forces acting on the rotor per unit width as 


where 



= B xl {y) 


= 5'-Ia>a.. 

„=a i 7i V *. 


(23) 


(24) 

(25) 


By — Byl(y) 






(26) 


The parameters in Table 1 were used in the following numerical examples. Figure 4 shows the flux distribution 
at zero speed and Figure 5 shows the convection effect and decrease in skin depth as velocity increases. Integrating 
the total flux within the air gap and within the journal, Figure 6 shows how the flux is compressed out of the 
journal as eddy currents increase. Fig. 7 shows the variation in lift and drag forces. Observe the change over to 
repulsion at 50,000 rpm. The drag force peaks at a value of 1100 Nt at fi = 10 8 rpm. 

The effect on the values of the forces of the number of terms used in the Fourier series was investigated. The 
magnetic forces were calculated at different spinning speeds using 10, 25, 50, 100, and 500 terms in the series. Figure 
8 shows the variation of the lift force. There is very little change in the results when the number of terms used in 
the series is 50 or more. There is approximately 3% change in the lift force when the number of terms is increased 
from 10 to 25 or from 25 to 50. This large number of harmonics is reasonable if the Fourier approximation to the 
applied square wave is plotted. The drag force effectively converged within 10 harmonics. The authors typically 
use 25 harmonics in dynamic analysis. 

An alternative approach is developed in [2] to find net magnetic field for all four magnets as a single system 
(simultaneously). The problem is posed in Fig. 9, which shows all four B fields. This complex boundary condition 
is expanded in a Fourier series and the previous equations applied. However, this approach requires the assumption 
that all the gaps are equal. At low spinning speeds there is no difference between the two methods as shown in Fig. 
10. Only at very high speeds do the two differ. (For parameters in this paper, 4% at 10 5 rpm). Researchers may 
wish to pursue this method. 



CONCLUSION 

A very general analytical expression for bearing forces has been presented. The model is applicable to journals 
made of permeable or nonpermeable material and to DC fields or AC fields. Both kinds of eddy current repulsion 
bearings are included. The predicted loss of lift (about 30% at 12500 rpm), agrees approximately with experimental 
results of Yamamura [11]. The expression for the drag force under a single sinusoidal field density wave is the same 
as that obtained by Meisenholder [9], 

Without considering eddy current effects, there is no speed dependent term in the constitutive model for a 
magnetic bearing. Eddy currents cause a loss of effective lift which can be viewed as an external load and may 
cause a change in equilibrium point, and change bearing coefficients. More importantly, eddy current drag causes 
coupling between (ar, y). Any change in the x gap affects the drag force in the y direction. Other work has shown 
this coupling effects bearing stability. 

A final comment applies to bearing coefficients (equivalent stiffness and damping matrices). The situation 
for a magnetic bearing is different from that of journal bearings. It is not possible to determine the 8 classical 
coefficients just by differentiating the force expressions with respect to £,^,£,77. The design of the control system 
may introduce state variables within the controller. In that case the dynamics of the electrical components must also 
be incorporated. Some subspace reduction might be possible, but it is impossible to determine even approximate 
coefficients without a completely designed control system. 
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NOMENCLATURE 

a Distance of magnet pole corner (near) from center line 

A Area of magnet pole face. 

b Distance of magnet pole corner (far) from center 

B Magnetic flux density 

Sq Magnetic flux density under the magnet pole 

D Diameter of journal 

e Unit vectors 

Fmag Magnetic force 

H Magnetic field intensity 

j V-t 

k Wave number 

L Inductance of circuit 
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Li Flux path length 

ni,n 2 Real and imaginary parts of 7 

N Number of turns in coil 

q Complex wave number 

R Resistance of coil 

Si Nondimensional velocity 

t Time 

V Velocity vector 

V Velocity (magnitude) 

Wd Field energy density 

Wj Field energy stored 

£ shaft displacement 

rj shaft displacement 

7 Complex number = cosh kA + Hoq/fi^ sinh k A 

A Gap length 

T Magnetomotive force 

6 Magnet pole angle with respect to magnet axis 

A Flux linkage 

fly fi r Permeability 

a Electrical conductivity of journal material 

E Closed surface, Summation sign 

$ Magnetic flux 

Q Spin speed of rotor 


APPENDIX A 

ESTIMATION OF MAGNETIC LIFT FORCE USING CIRCUIT THEORY 


If eddy currents are neglected, magnetic circuit theory can be used to approximate the magnetic lift force of 
a single magnet. The following assumptions are made in deriving the expressions for the magnetic lift force: 

1. Field fringing is neglected. 

2. Magnetization curve is linear ( B = fiH). 

3. Magnetic flux density B and field intensity H are uniform over cross-sections of the core, gap, or 
journal. 

An electromagnetic circuit is considered whose elements are gap, core, and journal. Each element has constant 
crossection Ai and length Li. The magnetic flux <j> is assumed constant throughout the circuit, and T is the total 
magnetomotive force within the circuit elements. 

The density relationship <f> = BA and the constitutive law B — fiH can be used to express field intensity in 
the ferromagnetic material in terms of the field intensity within the air gap. 

By definition T = / Hdl Integrating around the circuit and equating T to the current linked (NI) results in an 
equation for field intensity within the air gap. The field energy is determined within each element w = ^fiff and 
the total field energy is obtained by a volume integral over all the elements. 

By definition, the force is the rate of change of stored field energy with respect to the mechanical displacement. 


Lmag — 


fi 0 N 2 AI 2 


(2x + + ^-) 2 

V Prl Jir3 } 


(A.l) 


In addition, the total flux </>, the magnetic field density B, the magnetic flux linkage A, and the inductance L can 
be expressed in terms of the gap length x and current I as: 


4> = 

B = 
A = 


fioNIA 


(2x + i*- -f "-) 

Ho NI 

v T /in T Mrs' 

fi 0 N 2 IA 

(2x + 

' Vrl ^ l*r2' 


(A-2) 

(A.3) 

(A.4) 
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TABLE 1 

Four Magnet Bearing Parameters 


Diameter of journal D 
Path length in journal L 2 
Relative permeability of journal material fi v 2 
Electrical conductivity of journal material a 


Pole angle 

Path length in core L\ 

Relative permeability of core material ji r \ 

Area of pole face A 

Width of pole face W 

Number of turns N 

Distance of pole corner a 

Distance of pole corner b 

Resistance of coil R 

Permeability of free space /Jo 


15 degrees 
0.28 m 
10000.0 
0.00025 m 2 
0.008333 m 
800.0 
0.0075 m 
0.0375 m 
3.0 ohms 
Air E - 07 H/m 
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Figure 6. Flux distribution between air gap and journal. 
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